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ABSTRACT 


A formalism is developed for temper ature- dependent, self-consistent 
phonons in quantum solids with defects. Lattice vacancies and interstitials 
in solid helium and metallic hydrogen, as well as electronic excitations in 
solid helium, are treated as defects that modify properties of these systems. 
The information to be gained from the modified phonon spectrum is 
discus sed. 
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I. INTRODUCTION 


In this report we are concerned with phonon properties of quantum 
solids with defects. We shall mainly be interested in two substances: 

(1) electronically excited solid helium, and (2) metallic hydrogen. The class 
of defects that will be considered consists of electronic excitations in solid 
helium and lattice imperfections (vacancies and interstitials) in both solid 
helium and metallic hydrogen. The two substances considered here are of 
great current interest, both scientifically and because of potential practical 
applications. In particular, the long lifetime (~2. 3 h) and the high excitation 
energy (~20 eV) of the metastable 2 S state of atomic helium suggest the 
possibilities that high concentrations of these states might be achieved, under 
suitable circumstances, in the liquid and solid phases of helium, and that 
they could be utilized for storing energy (Ref. 1), The collisional quenching 
of the excited 2 S state is expected to be less severe in the solid than in the 
liquid, thereby singling out solid helium for primary consideration. 

The existence of metallic hydrogen was first proposed by Wigner and 
Huntington (Ref. 2). It is presently believed that the bulk of the interior of 
Jupiter is metallic hydrogen, most likely in the liquid phase. An added 
impetus for the study of metallic hydrogen is the expectation that it is a high- 
temperature (100 to 300 K) superconductor (Refs. 3 and 4). 

Quantum solids are characterized by the fact that the zero-point 
motion of constituent atoms occurs over a sizeable portion of the atomic 
cell. It is well established that solid helium is quantum. Although metallic 
hydrogen has also been assumed in the past to be quantum, recent work by 
Caron (Ref. 5) seems to indicate that this assumption is not entirely proper. 

It appears that metallic hydrogen falls somewhere in between classical and 
quantum solids; we shall refer to it as semiquantum. Special techniques 
have been developed to treat the highly anharmonic quantum solids, since 
standard approaches for classical solids, which are based on the harmonic 
approximation, simply break down for quantum solids, as manifested by 
phonon frequencies becoming imaginary (Ref. 6). Although classical 
methods may be marginally valid for metallic hydrogen, it is felt safer to 
treat it by quantum- solid techniques, as we shall do in this report. 
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As a result of the large zero-point motion of atoms in quantum solids, 
"c _ ertai'n — latti'c'e - d-ef'e'clrs^ - s'urc‘h _ a"s _ va"C'aiTci"e"s _ a‘n'd _ i'nt'e'rstiti'a"l"s^ a‘r _ e — e _ a"S'il'y 
excited, even at zero temperature* These defects are expected to have a 
marked effect on various physical properties of a quantum solid, especially 
on those that depend on the lattice structure* In particular, the phonon 
properties are expected to be strongly influenced by the defects. It is the 
primary purpose of this report to discuss phonons in quantum solids with 
defects, in the context of solid helium and metallic hydrogen. In the case of 
metallic hydrogen, the important defects we shall consider are vacancies 
and interstitials. In the case of solid helium, our primary interest is elec- 
tronic excitations, treated as lattice imperfections; however, our discussion 
will be sufficiently general to include also the '’normal' 1 lattice defects, 
vacancies and interstitials. 

Most of the work on defects in lattices has been done for classical 
solids (Ref. 7). Theoretical studies of defects in quantum solids were initi- 
ated around 1970, and have been mostly concerned with two problems. The 
first, in analogy to the classical case, is concerned with determing the mod- 
ifications to the phonon spectrum due to various types of defects. Varma 
(Ref* 8) considered a single isotopic substitution, and Aksenov (Ref. 9) 
studied finite concentrations of static vacancies* Our intent is to generalize 
the approaches of these papers to treat situations relevant to electronically 
excited helium and metallic hydrogen at finite temperatures. The second 
problem that has been studied in the past is concerned with using a 
collective-excitation approach to treat lattice defects and their motion in 
quantum solids (Refs. 10 through 13). In this report we shall not treat this 
problem. 

II. SELF-CONSISTENT PHONONS IN A PERFECT QUANTUM SOLID 

To find the phonon spectrum of a system with a not too high defect 
concentr ation it is very convenient to use the phonon spectrum of the corre- 
sponding perfect and undistorted crystal* The periodicity of the latter can 
be used advantageously to define normal coordinates and greatly reduce the 
diagonalization problem* We review here the basic results on self- 
consistent phonons (for exhaustive reviews, see Refs. 14 through 16) in a 
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form that will be convenient for treating the defect problem. To introduce 
the notation we first briefly summarize the derivation of phonons in a per- 
fect, classical system (Ref. 17). 

A. CLASSICAL LATTICE DYNAMICS 

The basic assumption of small vibrational amplitudes allows one 
to perform a series expansion of the potential in the displacements and to 
truncate the series after the dominant, quadratic terms. The original 
Hamiltonian is then replaced by the approximate, harmonic Hamiltonian 


wher e 



M. u 2 

i i a 


iat 


1 V * ic * 

+ o* 7 $ .q u. u.o 

2 # jP ux j(3 

ice 

jP 


_ ia I a 2 V 

" \ d *ia 


eq. 


(1) 


are the force constants derived from the original potential, and a, p = x,y,z. 
The classical equations of motion of this Hamiltonian are, of course, linear; 


M. u. 
i i a 


jp 


( 2 ) 


The index i (or j) counts the atoms and is shorthand for the pair mp (or nv), 
where the Roman letter denotes the unit cell and the Greek letter a site in 
the cell (allowing for more generality than a Bravais lattice). The system 
shown in Eqs. (2) can be greatly uncoupled by assuming the following space 
and time dependence of the displacements: 


i [ k ' R m H i " W(k)t ] 


(3) 
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where the wave vector k reflects the lattice periodicity. Note that masses 

— M do - not depend' on the celTlndex m. The System shown in 'Eqs.TZ) now 
P- 

becomes 


J(k)U ]iQ (k) = X) D Jp(k>u v p(k) 

vP 


(4) 


where 


. _ ik- {R -R ) 

D^a(k) = I (5) 

V P IMTT £-> ‘"'P 

v p V n 

the so-called dynamic matrix, is the spatial Fourier transform of the force 
constant matrix. Let r denote the number of sites per unit cell. Compar- 
ing Eqs. (2) and (4), we see that the explicit introduction of lattice periodic- 
ity has greatly simplified the problem: while Eqs. (2) are a system of 3rN 

coupled equations (N being the number of unit cells in the system, practi- 
cally infinity), Eqs. (4) consist of N systems (one for each k) of only 3r 
coupled equations each. To each k there correspond 3r roots, denoted by 
\ = 1 , • • • , 3r. The corresponding eigenvectors of the dynamic matrix, 

e (k\) , the so-called polarization vectors, fulfill the orthogonality and 
pa 

closure relations 



Once the dynamic matrix has been diagonalized (which is easily done 
because it is low dimensional), the harmonic Hamiltonian can be shown to 
be completely diagonal in the normal coordinates Q(kX.), related to the 
atomic displacements by 
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Indeed, substitution in Expansion (1) gives 


( 7 ) 


H h = 1E 

kx 

which is, as expected, the Hamiltonian of uncoupled harmonic oscillators. 
This form is extremely important in the treatment of a quantum solid, since 
it allows one to write down explicitly approximate wave functions of the solid. 

B. VARIATIONAL PRINCIPLE 

The basic idea of self-consistent phonons in quantum solids, for which 
the atomic displacements are too large and hence the harmonic part of the 
Expansion (1) breaks down even as a zeroth- order approximation for a pos- 
sible perturbation treatment, has been to use a variational approach. 

Instead of expanding, one is looking for the most general harmonic 
Hamiltonian (Expansion 1) that would best approximate the given system in 
some sense. The most common approach, systematically and successfully 
developed by Koehler (Refs. 18 and 19), has been to regard the effective 
force constants, or equivalent quantities, as variable parameters determined 
by minimizing the ground-state energy with respect to them. The energy is 
evaluated by averaging the given, anharmonic Hamiltonian over harmonic 
oscillator wave functions. Consequent minimization gives a self-consistent 
condition for the phonons, with the frequencies appearing not only explicitly 
as in Eq. (4), but also in the dynamic matrix. 


Q(k\) 


2 + to(k\) 2 


Q(k\) 


(8) 
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Because we are interested in temperature-dependent self-consistent 
phonons (es pe cradly for metallic hydrogen), we have to^ replace the ground- — 
state energy by the corresponding thermodynamic quantity, the free energy 
F, given by 


e -PF . £ e -? E „ = Tr 
n 

where the summation is over all states of the system and |3 = 1/kT. The 
statistical operator, or density matrix, is given by 

e - PH 


The thermodynamic average of an operator O is 


<Q> 


Tr (pO) 


Tr (O e"^ H ) 
Tr (e- pH > 


The extension of the zero-temperature , ground- state variational principle 
to finite temperatures is the thermodynamic Gibbs - Bogolyubov principle 
(for a proof and review see, for example. Ref, 20), which states that 


F £ F trial = F 0 + < H - H n> 


0 


Tr 


s F 0 + 


-PH ' 

(H - H q ) e 


Tr 


•|3H 


0 


(9) 


where Hq is a zeroth order Hamiltonian, conveniently chosen as a 
Hamiltonian for which the problem can be solved explicitly. For T = 0 the 
first and last terms in Eq. (9) cancel and the principle reduces to the 
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well-known variational principle. This variational principle will be used 
now to derive temperature-dependent, self-consistent phonons. 

C. SELF-CONSISTENT PHONONS 

To obtain the self-consistent phonons we first evaluate the trial free 
energy in Eq. (9) and then minimize it. The Hamiltonian is of the form 


H 


T + V 


t 2 i a 2 

2 a 2 

mjiO' n mjj.O' 
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4 v ( ir + u 

2 / j I m(i mp. 


m[Ji 

nv 


- R 


nv 



( 10 ) 


and the zeroth order Hamiltonian, H q , is naturally chosen to be the har- 
monic Hamiltonian in Eq. (1), which has the advantage that all thermo- 
dynamic qualities and averages can be evaluated explicitly. The kinetic 
energy terms cancel, and we are left with 

F trial - F 0 + < V - V 


It should be noted that in the case of solid helium the potential V appearing 
in Eq. (10) is not the usual interatomic potential (of the Lennard- Jones 
variety), but an effective potential that incorporates short-range correlation 
effects (Ref. 21). The point is that the potential that a pair of helium atoms 
in the lattice see is drastically modified by the surrounding atoms, so that 
the pair potential derived from gas data is a very poor description of the 
interatomic forces in the solid. In the case of metallic hydrogen, V is the 
screened ion-ion potential (Refs. 22, 5). 

In terms of the normal coordinates (Eq. 7), which diagonalize the 
harmonic Hamiltonian (Eq. 8), the free energy is additive and the density 
matrix is multiplicative, and all necessary quantities are explicitly avail- 
able for harmonic oscillators (see, for example. Ref. 23 or Ref. 24): 
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The averages of V and are given by 



The result of Eq. (14) is obtained by averaging the square of the displace- 
ment of a single harmonic oscillator (see Eqs. 8 and 12). In the low- 
temperature limit each term gives "hco/4, which is the potential energy 
contribution to the zez’o-point energy. In the high- temper ature limit, each 
term contributes kT/2, the expected classical result. 

The trial free energy is given by combining Eqs. (11), (13), and (14). 
The resulting expression depends on all theto^^, which are obviously 
related to the effective force constants in Eq. (1). It is convenient to regard 
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these frequencies as variational parameters, to be determined by imposing 
the condition 


9F, . . 
trial 

9 w k\ 


0 


(15) 


This procedure is similar to that employed by Morley and Kliewer (Ref. 25) 
to derive zer o-temperatur e , self-consistent phonons. The algebra is some- 
what long, but straightforward, and the result is 



This is the basic result of the variational self-consistent phonon 
approach, replacing second derivatives of the potential at equilibrium by 
thermodynamic averages of the second derivatives over harmonic density 
matrices. This system has to be solved self- consistently , as the thermo- 
dynamic average depends on all the frequencies. To relate this result to 
the corresponding classical result, Eqs. (4) and (5), we transform from 
normal coordinates to atomic displacements by means of Eq. (7). The 
result is formally similar to Eqs. (4) and (5) except that the dynamic 
matrix now becomes 




/M M 
V V- v 



ik* (R 

nv 


- R ) 

mp 


(17) 
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Both this and the classical dynamic matrix can be written in more 
— — —detail —specifying— prope r-tie s-of acoustic-and optic modes,— by -writing-the. 
potential as a sum of pair potentials, as follows: 


7 

d V 


a 2 
9 v 

mp, nv 

-an * nv 

m jjl at nvp 


9 U 9 U o 

mpcr nvp 


*2 
O V 


m p, rp 


rp 


9u 9u a 

mpa' mpp 


r np - nv 


Substituting this into the dynamic matrix, we get 


( 18 ) 



(19) 

D. EVALUATION OF THERMODYNAMIC AVERAGE 

The last step is to evaluate the thermodynamic average of the second 
derivative of the pair potential, and to reduce the multidimensional integral 
to a form that can be calculated. This is achieved by Koehler 1 s method 
(Refs. 18 and 19) of expressing the density matrix Eq. (12) in terms of 
atomic displacements (resulting in correlated Gaussians) and consequent 
reduction of the multidimensional integral appearing in the thermodynamic 
average to a three-dimensional one, which can be evaluated numerically. 

We start by substituting Eq. (7) in Eq. (12), and Eq. (12) in Eq. (19), with 
the result 


*2 
9 v 

mp, n v 

9U 9U ry 

mp mppy 



9 2 v 


mp, nv 


-£ 

rpy 

s<r6 


G r ‘ u u 
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m par mpp 


( 20 ) 
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where the temperature-dependent correlation function is 


G rp I 

s<r o 


N 


E N /MpM < 
kX 


^kx 


kx ^ 

tgh 2kT 


(KX) 


CTO 


(21) 

Next, since the second derivative depends on the coordinates of two 
particles only, we can perform the integration over all coordinates but these 
two. For this it is convenient to partition the matrix G and its inverse 
G ^ into a six-dimensional block involving the indices of particles m p, nv, 
and the rest: 


mp nv 



The result of the integration is (Ref. 18) 


8 2 v 

mp, n v 

~9u 9u p 
mpa m[i.p y 


(det h l) 


7 T 



a 2 v 

m p, nv 
mpQ' mpp 


Uh" 1 U 


where U is a six-dimensional vector consisting of the components of u 

_ y _ ^ nap 

and u . Since G and G are symmetric, h and h have the block 
nv 

structure 


with 


h. Vi h-‘.(v 

\. B i A/ \ D 1 C 


C = (A - BA -1 B T ) , D = (B T - AB -1 a) 
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For central force potentials the integral can be further reduced by 
introducing relative - and center of ma^ 2 “oo'r‘dinates 


U 


m \i 



u 


nv 



u h" 1 u = 2R(C + D)R+-|r(C-D)r 


The integration over R is readily performed: 




a 2 v(|a 


mp 


9r 


a 





r{C - D)r 


( 22 ) 


This form is useful if the three-dimensional matrix (C - D) can be related 
to the correlation matrix G (Eq. 21). Some matrix manipulation shows that 
C - D = (A - B)” 1 . A and B are the diagonal and off-diagonal blocks of the 
inverse of the correlation matrix G. Fortunately, the latter can be inverted 
(despite its practically infinite dimension) by using the orthogonality and 
closure properties of the polarization vectors. It is readily verified that 



m|ia 

nup 


-Y' 

N Z-/ 


iiaj 




k\ 




ctgh 


IkT V <kM {kM 6 


ik- (R - R ) 

nv m|j. 


(23) 


Defining the matrix F = (C - D) ^ = A - B with elements F(mp, nv)^ 

= (G - (G and collecting all the results, we obtain the self- 

x mpp nvp * 

consistent phonon scheme 


“L v (kX| ■ £ 

r P 

vp 



12 


JPL Technical Memorandum 33-701 





/ ik ■ (R -R )\ 

x (e^kX) - e p(3 (k\) e r P m Pj (26) 

This is the self-consistent harmonic approximation in a form equiva- 
lent to that of Gillis (Ref. 26), and will be used as a basis for treating the 
defect problem. 


III. SELF-CONSISTENT PHONONS IN AN 
IMPERFECT QUANTUM SOLID 


We formulate a scheme for computing self-consistent phonons in 
quantum solids with defects in the expectation of applying it to systems such 
as solid helium and metallic hydrogen, in which electronic excitations and 
lattice imperfections are regarded as defects. The first system is of 
interest in connection with the possibility of energy storage in metastable 
excited states (Ref. 1). The phonon spectrum for metallic hydrogen was 
evaluated (Ref. 5) for a perfect crystal at zero temperature. It is important 
to obtain the temperature-dependent phonon spectrum, including vacancies 
and interstitials that set in with temperature, and to see how the defects 
modify the phonon spectrum and affect the lattice stability and, possibly, 
melting of metallic hydrogen. Moreover, lattice defects are expected to 
influence the superconducting properties of metallic hydrogen, as will be 
discussed in Section V. For solid helium the problem is similar to the pas- 
sage from perfect to imperfect classical solids (for an exhaustive summary 
see Ref. 7), only that now a quantum solid is being considered. The two 
types of defects, electronic excitations and lattice imperfections, differ 
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greatly in their activation energy (20 eV and fraction of an eV, respectively) 
and hence deserve to. i>_e -tne-at&d- d-iff-e-r-ently; _ ' The fact that lattice imperfec- 
tions have low activation energies means that their concentration is greatly 
temperature-dependent and so requires a statistical-mechanical treatment. 
This is not so in the case of electronic excitations, which cannot be ther- 

5 

mally excited because their effective temperatures are of the order of 10 K, 
enormous compared to the solid-helium temperatures of £ 10 K. For this 
reason we shall treat electronic excitations as fixed, temperature- 
independent impurities with a random distribution in the solid, characterized 
only by an average concentration. 

A. VARIATIONAL PRINCIPLE 

The variational principle for the free energy (Eq. 9) was stated for the 
canonical ensemble, in which the number of particles was fixed. To allow 
for a variable number of defects we have to use the grand-canonical ensemble, 
in which the number of particles (defects in our case) is allowed to vary. 

The statistical quantities we need are defined as follows (Ref. 27). The 
grand partition function is 


Z 


G 


Tr e 


■P(H - p N 
'v v 


W 


(27) 


where Tr means summation over the states as well as over the number of 
vacancies and interstitials, and \x and |jc are the chemical potentials asso- 
ciated with these defects* The analog to the free energy is the thermo- 
dynamic potential Q, 


£2 = -kT in Z, 


while the density matrix, or statistical operator, is given by 


1 -KH-„ v N v - h .N.) 
P G = ~zZ e 


- e 


p(f2-H + p v N v + p.N.) 


(28) 
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The average number of vacancies, for example, is given by 


<N v > = Tr (N v p G ) = 


Tr N e 
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-P(H - u N - u.N. )' 
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( 29 ) 


and, similarly 


an 

9 ^, 


<N.> 


an 

ap.. 


(30) 


A variational principle will now be derived by relating the grand 
canonical to the canonical quantities and using the variational principle for 
the latter. For a system of N particles we have (Eq. 9): 


F < TT 

N “ N, trial 
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The grand partition function is related to the canonical partition function 

Z N by 
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Using the variational upper bound on we have 


-pS2 , C"> P(F N, trial " "^trial 

- 2^ e 
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This defines the variational principle for the thermodynamic potential: 


o o 1 , Y' " P(F N, trial " ^ 

52 " ^trial ~ ~B fn / , e 
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= 2 e 
N 


-P 
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(31) 


In the case of our defect problem, N is the set of numbers of vacancies and 

interstitials, N = <N , N. 

( v 1 

will be differentiated with respect to phonon frequencies to obtain self- 
consistent phonons, with respect to effective equilibrium positions to find 
the lattice relaxation due to defects, and with respect to the chemical poten- 
tials to find the defect concentration (Eq. 29). 


and uN = \x N + 
^ r v v 


p.N.. 
r i l 


In what follows 


trial 


B. DEFECT HAMILTONIAN 

The first step in evaluating the thermodynamic potential is to write 
down the Hamiltonian of the system in the presence of defects. The 
Hamiltonian of a perfect crystal is 


H 


i ij 


(i ^ j in the summation) 


In the presence of the defects enumerated above (electronic excitations 
and/or vacancies and interstitials) the Hamiltonian will be modified in two 
important ways. First, the presence of electronic excitations drastically 
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changes the interatomic potential between two atoms, one of which at least 
is excited* The excited sites will be characterized by 



! 1 site i excited 
0 site i normal 


The modified potential becomes 


u. . 


v. . + 6v. . 
*3 1J 


(,<*> 


(x; 

<T 

J 




+ 6 V 




J 


(32) 


where 


6 v. . = vP - v.. 
ij ij ij 


c, XX 

O 1 V. * = V, . - V. . 

13 1J 13 

v gx is the potential between two atoms, one of which is electronically 
xx 

excited, and v is the potential between two electronically excited atoms* 

CT XX 

For helium, for example, both v & and v" are available from extensive 

quantum- chemical calculations (Refs* 28 and 29)* The potential u„ was 

constructed in such a manner that when one atom is excited u. , = v§ x , and 

11 

xx * 

when both are excited u. . = v 

*3 

To describe the effect of missing or superfluous atoms we similarly 
introduce 



1 i occupied 
0 i vacant 



1 i 1 vacant 
0 i T occupied 
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, (v) = 

( 0 i 

occupied 

i 

( 1 i 

vacant 

(i) _ 

(0 i' 

vacant 

i 1 

1 1 V 

occupied 


The superscripts (v) and (i) stand for vacancies and interstitials, and the 
primed index runs over interstitial sites. The modified Hamiltonian can be 
written as 






1 V (i) (i) ^ 1 V* 

rr / u. cr;. <r . + -=■ / U. . , 

2 Z — 4 i 1 J ’ i' j' 2 lj ' 


.(v) 


+ 2 


zE 




ij 1 


p< V > 

1J 1 J 


I'J 


The first two terms are the kinetic energies of the occupied regular sites 
and of the interstitial sites, respectively. The third term is the potential 
energy between the regular sites and the fourth term between two inter- 
stitials. The last two terms describe the potential between one regular 
site and an interstitial, and are obviously equal to each other. This way 
of writing the Hamiltonian is very similar to that of Aksenov (Ref. 9), who 
treats vacancies in quantum solids. Our work differs from his in three 
important aspects. First, he considers vacancies only, and no interstitials. 
Secondly, he uses the canonical distribution throughout, while it seems to 
us essential to use the grand- canonic ensemble when the number of particles 
varies. Third, Aksenov does not include lattice relaxation around defects 
that can be incorporated in our formalism. Such lattice relaxation effects 
are expected to be very important in solid helium. 

We next express the p's in the Hamiltonian in terms of the correspond- 
ing cr's and partition the Hamiltonian into lattice and defect parts: 


H 


H 


latt 


+ H 


def 


(33) 
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( 34 ) 


H latt " 2 *i + 2 ]C U ij 

i ij 


H 


def 





+ 



a. . 


? 


(■ 



(V) (v)\ 

°i / 


(35) 

The defect Hamiltonian is the or-dependent part of the Hamiltonian and 
obviously vanishes for c r = 0 3 i. e. , in the absence of defects. 

In the variational principle (Eq. 31) we choose as the harmonic 
part of the lattice Hamiltonian, so that 


<H - H > = <H 

0 


latt - H 0> + < H def> 


(36) 


The "perturbation” Hamiltonian consists of two parts, the effective anhar- 

monicity of the perfect crystal and the average defect Hamiltonian. As we 

will see, the first term contains the self-consistent phonons previously 

derived, as well as the modification due to electronic excitations (through 

the modification of the potential, Eq. 32). The second term contains the 

effect of vacancies and interstitials. The average in^H, consists of two 

def 

parts, the average of the kinetic and potential energy operators over the 

density matrix of harmonic oscillators, and the summations of the linear 

and quadratic terms in the tr ' s over N and N. . As a result we have 

v 1 
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. < H det> = 



- 0 * 


(V) 


)_<t> + jr 



iC <Vj'> 

i'j' 


_(i) 


(v). 

(T ) 


£ <v - i ^ 

ij’ 


(2 - 




<i> 


(37) 




i (v) , (i) 

where cr and c r 


are the average vacancy and interstitial concentrations. 


It should be pointed out that in passing from H^ e £ to ^H ( £ e £^> one 
obtains a translationally invariant defect Hamiltonian, thereby suppressing 
all information about the localized nature of the defect, A more rigorous 
description of a finite defect concentration would be provided by a scheme 
such as the coherent potential approximation. The discussion in Subsection 
C below is kept general in that translational invariance of the lattice is not 
assumed. Calculations will be carried out in the low-defect concentration 
approximation when only linear terms in (r( v ) and cr(i) are retained 


< H def> “ 


(Or 


:i) (V 


ff ( V )) <T>+ ,(i)^ <Uij| >. 

ij’ 


(v) 


<u. > 

N ir 


ij 


N. - N 

l v 

N 


<T> 



ij 1 



ij 
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This form is next substituted in the variational principle, Eq . (31), resulting 


in 


^trial 




N , N. 
v* 1 


-« F o * - v 






(39) 



This derivation seems somewhat inconsistent. On the one hand, only linear 

terms in N and N. were retained due to the assumption of few defects. On 
vi 

the other hand, the summation was performed over all values, from 0 toco, 
using the well-known sum of a geometric progression. In the Appendix we 
start from the total defect Hamiltonian Eq. (37) rather than Eq. (38), derive 
an integral representation for the trial thermodynamic potential, and show 
that in cases of interest to us the result reduces to Eq. (38), 


The first two terms in Eq, (38) describe the variational free energy of 
a crystal with some fixed, localized electronic excitation, and will be con- 
sidered in the next subsection. The last two terms are the contribution of 
interstitials and vacancies, and will be considered later. 


C. ELECTRONIC EXCITATIONS AS DEFECTS 
The first two terms in Eq . (39), 

F trial = F 0 + < U - V 
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are of the same form as the variational free energy of a perfect crystal. 
Denoting the normal modes by , we— obtain from the previous- section, by — 
minimizing the free energy with respect to the normal mode frequencies, 
the basic equation of lattice dynamics 



Despite the formal similarity, Eq. (40) differs from the equation for a per- 
fect crystal in that the potential U is modified close to electronically excited 
sites, as described by Eq. (32), As a result of different potentials between 
different pairs of atoms, there is lattice relaxation around the defects. The 
lattice loses its periodicity, and we can no longer use the translational 
symmetry to simplify Eq. (40). The normal modes differ from those 
of the perfect crystal, both in their values and in the appearance of localized 
frequencies not present in the perfect crystal. 

The most convenient approach to Eq. (40) is to write U as the perfect 
crystal potential V plus a correction term, thus enabling one to use the 
well-known solution of the perfect crystal as a basis for the defect problem 
(Ref. 7). As we shall shortly see, the correction terms in the resulting 
equation for the normal modes can be characterized by a defect matrix, 
consisting of deviations in mass and in force constants (i. e. , interatomic 
forces) from the corresponding values of a perfect crystal. For a quantum 
solid the first type of defect has been discussed by Varma (Ref. 8), who 
considered a single isotopic substitution, while we consider here a change 
in the interatomic potential. Although we expect more drastic effects due 
to electronic excitations (not unlike bubbles in liquid helium with electronic 
excitations) than due to an isotopic substitution, the formal approach to both 
problems is rather similar, as we now proceed to show. 

Looking at electronic excitations, we can consider either a single 
excitation at the origin or a finite concentration of randomly distributed 
electronic excitations. In the first case the formulation becomes practically 
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identical to that of Varma (Ref. 8). Of considerably greater practical inter 
est is the second case, which we now examine. By Eq. (32) we can rewrite 
Eq. (40) in the form 


Mw! 



= 2cr ( 1 - a - ) 


a z 5v 

\3q v 3c h 


+ (T 


a 2 v 
\9q v 9q v 


(41) 


where a is the concentration of electronic excitations* The left-hand side 

looks like the expression appearing in the perfect crystal equation, but this 

is not quite the case. The potential is evaluated at relaxed lattice positions, 

R + 6R + u, where the R are the equilibrium positions of atoms in the 
o o 

undistorted lattice, 6R are the lattice relaxation displacements due to the 
defects considered, and u are the phonon displacements. We can write 


V(R n + 6R + u) = V(R q + u) + V(R n + 6R + u 


0 


V(R Q + u) 


= V(R q iu)+AV = V nd + AV 


where AV is the distortion potential due to lattice relaxation, and is 

the undistorted lattice potential. Equation (41) can now be written in the 
form 


Mcj 


Z 

v 



9 Z AV 

v 9 q v 9q, 


+ 2 cr ( 1 - o ) 


8 2 5V 

v 9 q^ dc i\ 


+ a - 


8 Z 6'V 

, 9 qy 8 q, 


(42) 


It is convenient to transform from normal modes to atomic 
displacements 


u, 

ia 



v 


T. 
ia v 


q 


v 
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Because of the lattice distortion the do not possess translational 

--- . symmetry: an.d_.we_c.an not introduce p olaxizLati on vector s , as we did for— per- 
fect crystals. The have the orthogonality and closure properties 



6 


v v 1 



T T 
iav j(3v 


6 .. 

U 



The second derivatives in Eq. (42) become 


d = V T T* a - 

itj|3 ^ iP,,3u - aU jP 


Using the orthogonality, Eq, (42) can be rearranged to read 



where C is the defect matrix for the problem 


C 


ic* 

jP = 



'6' V 



(43) 


(44) 


For a vanishing defect matrix the problem reduces to that of a perfect 
crystal. The defect matrix consists of two parts, one due to the modified 
interatomic potentials, and the second due to the deviation of the distorted 
lattice potential from the undistorted one. 
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Equation (43) is conveniently solved in terms of classical Green 
functions, defined as solutions to Eq, (43) with the right-hand side replaced 
by a 6 -inhomogeneity, 


• Mu 


IQ 

Slv 


(w v ) + 



(“v) 


6 


io? 

iv 


It can be seen, by direct substitution, that the Green functions are given by 


g jP 


(“v) 

k\ 


V < kx > T ,|i < k ^> 



2 


V 


(45) 


where are the normal modes of the perfect crystal, and T^(k\) are the 

corresponding eigenvectors, related to the polarization vectors by 


t. (k\) 

lot 


* 

Am ° 


(k M e' ik ' R 


Knowing the Green functions, it can again be verified by direct substi- 
tution that the solution to Eq. (43) can be written as 


T iov = £ <“v) CJP T lyv 

jP 

w 

This can be regarded as the matrix equation 

( 1 - g C) T = 0 

solved by diagonalizing the product of the Green and defect matrices 


det 



gj(3 <“>v) 



0 


(46) 
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Unlike the secular equation of a perfect crystal, whose size is determined 

geometric_ally_by the number of vibrational degr_ees of free.dormper unit cell, 

the size of this equation is determined physically by the number of neighbors 
affected by the defect* 

Finally, the lattice distortion 6R can be determined by varying the 
free energy with respect to 6R: 

F trial = F 0 - < H 0> 0 + < H latt < R 0 + 6R + “> + H def < R 0 + 6R + u '> 0 


Only the last two terms include the lattice distortion, and imposing the 
condition = 0 results in an equation similar to that derived by 

Varma (Ref, 8), 


D* VACANCIES AND INTERSTITIALS 

For this problem we have to include the last two terms in Eq. (38) 
that were not considered so far* The self-consistent phonons will be deter- 
mined from the condition 


B Srial 

Sw v 


0 


We conveniently define 


E. = 

l 
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Differentiation gives 


3ft 


'trial 


3F 


trial 





J_ / 9<T> 8 

N I 3oj v 8u> v / 4 

\ ij' 




N 




The derivatives required here are all available in the previous section for 
the perfect crystal. They are 


trial _ 

4 



(3hu) 

v 

2 


sinh 


^ v 
2 



9<T> 

du) w 



Ptiw 

V 

2 


sinh 


2 

2 




U. ^ denotes the interaction between regular lattice sites and interstitial 
sites. Defining the function sinhc x ~ (sinh x)/x, we collect and rearrange 
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all terms. The resulting self-consistent phonon equation in the presence of 
defects is __ . 


a 2 u \ z_ 

+ N 




E i /^int 


a 2 u 


1 - e. 


i Vq v H 


1 - E 


9( V 


rriGJ,, = 


, sinhc ptkw - 1 / E. 

i4 11 1 1 


(47) 


v 


N sinhc ptiio,, + 1 \ 1 - E, 1-E 
v \ l v 


We next evaluate the average number of defects from relations Eqs. (29) 
and (30): 


3 fb . , E 
<n > = - — iv al - v 

v 3u 


V 


1 - E 

v 


(48) 


9^2 . . E. 

<V= 

'l 1 


(49) 


Substituting these results into Eq. (47), we can write the self-consistent 
phonon equation in terms of the vacancy and interstitial concentrations: 




1 + (<r 


(i) 


, v sinhc Shu, - 1 

<r< v )) _ _ 

' sinhc (3 hco v + 1 



sinhc phw 
sinhc phw 


v 


v 


- 1 

TT 


\ 


(50) 
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This equation can be treated by the classical Green function technique 
described in the previous subsection. The quantities appearing here, par- 
ticularly the defect concentration, depend on the chemical potentials and 
P v that are not yet determined. One way to determine then is to compare 
Eq. (48) and Eq. (49) with 

-AE /kT -AE./kT 

<N v > = e V <N.> = e 1 

where AE^, AE. are the activation energies of vacancies and interstitials, 
and thus relate the chemical potentials to the activation energies. 

IV. SOLID HELIUM 

In this section we discuss the relevance of phonons to the study of 

electronic excitations in solid helium. Our attention will be exclusively 
3 

devoted to the 2 S metastable excitations of helium because only they have 

a sufficiently long lifetime to be studied conveniently experimentally. In 

fact, as noted in the Introduction, the long lifetime and the high excitation 

3 

energy of the He(2 S) s He state provide the major motivation for the study 
of electronically excited helium, for eventual energy- storage applications. 

A. -MEASUREMENT OF He* CONCENTRATIONS 

An excited atom in the solid-helium lattice can be regarded as a defect 
or an impurity. As such, it produces a static lattice distortion in its imme- 
diate neighborhood and thereby leads to a modification of the phonon spec- 
trum. The lattice distortion is expected to be quite substantial, in analogy 
to the case of liquid helium where M bubbles n or cavities around excited 

o 

atoms are known to exist, with estimated bubble radii of the order of 7 A, 

Q 

as compared to the mean interatomic distances of about 3. 5 A (Ref. 30), 

Both the static lattice distortions and the vibration of normal and excited 
helium atoms about the distorted equilibrium positions are treated self- 
consistently , at any temperature, by the formalism of the last section. 
Knowledge of the phonon spectrum of an imperfect solid-helium crystal, 
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together with sound attenuation measurements in the crystal, may provide 

* 

a viable technique for determining the concentration of He_ atom s_ in solid _ 
helium, as we shall now discuss. 

It is known that the presence of impurities in an otherwise perfect 
crystal leads to scattering and attenuation of sound waves propagating in 
the lattice. For sufficiently low-impurity concentrations, attenuation of 
sound waves is essentially proportional to the impurity concentration, 
although nonlinearities are expected to complicate matters at higher con- 
centrations. Analogously, by measuring the attenuation characteristics of 
sound waves in an electronically excited solid-helium crystal, one is 
expected to obtain information about the excited- state concentration. Exper- 
imental techniques for making such measurements on ordinary crystals with 
impurities are well established. There appear no major obstacles in apply- 
ing them to solid helium. To interpret sound- attenuation experiments, one 
must know the effect of excited- state impurities on the propagation of sound 
waves in solid helium. The work on the phonon properties of solid helium, 
reported here, is being extended to deal with this problem. 

Present experimental techniques for measuring sound attenuation in 
imperfect crystals are quite accurate and the errors in determining per- 
centage impurity concentrations are of the order of 0.01, for concentrations 
not exceeding 10%. It is expected that comparable accuracies can be 
achieved in the solid-helium case. It should be stressed that at present 

there is no reliable method of estimating excited- state concentrations in 
solid helium (Ref. 31) and that therefore sound- attenuation techniques are 
expected to find an important application in this area. 

B. DEEXCITATION PROCESSES 

According to theoretical estimates (Refs. 32 and 33), an isolated He 
atom has a lifetime of some 2.3.h. In a solid-helium crystal, the lifetime 
of the He atom is generally expected to be shorter in view of various per- 
turbations acting on the atom that tend to facilitate its transitions to the 
ground state. A case where an analysis of lattice perturbation effects has 
been performed is that of the excited N = N(2p D) state in a lattice of 
molecules (Ref. 34). The results of this analysis show that indeed a rather 
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drastic reduction of the N lifetime is caused by lattice perturbations. The 
case of solid helium is sufficiently different to require independent consider- 
ations, although one might a priori expect lifetime shortening effects here, 
too. 

The mechanisms by means of which an He atom may lose its excita- 
tion are quite numerous,. Some of the more important ones are: 

(1) increased spin-orbit mixing (and thereby enhanced radiative decay) 
caused by the presence of neighboring atoms, (Z) He migration in the 
crystal and subsequent M collisional ir He -He deexcitation, (3) formation of 
He -He molecules or other complexes having shorter lifetimes, (4) radia- 
tionless deexcitation into phonons, and (5) loss of electronic energy in 
exciting vacancies and interstitials in solid helium. 

To be able to discuss these various processes quantitatively, one 
must know how electronic excitations modify the static solid-helium lattice 
and how they couple to the lattice vibrations, i. e. , phonons. The formal- 
ism presented in Section III has been developed for the purpose of supplying 
this information. Here we shall briefly outline how the static and dynamic 
properties of the lattice enter into the calculations of some of the above 
decay processes. 

The fact that an He atom has spin 1 and that, apart from spin- orbit 
terms, the Hamiltonian is spin independent means that the only way He in 
a solid-helium lattice can decay is either by the action of spin-orbit forces 
or by mutual deexcitation with another He* atom. Spin-orbit perturbations 
lead to both radiative and radiationless decays. In the case of radiationless 
decays the electronic excitation energy of an He* atom is dumped into 
phonons to appear ultimately as heat. For energy- storage applications, it 
is extremely important that such radiationless He decays have very low 
probability. Estimates of matrix elements for these decays, at various 
temperatures, will be made on the basis of numerical calculations of self- 
consistent phonons. 

Although radiative decays seem to represent another clear-cut loss 
mechanism of He atoms, it must be borne in mind that phonons emitted by 
the decay of He atoms may become trapped or reabsorbed as they propa- 
gate through the crystal, so that, effectively, energy is retained in a useful 
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form. In fact, the photon concept may lose its meaning in a crystal, since 
photons are known to couple to electronic excitations (or excitons) in crys- 
tals to generate new types of excitations, the polaritons (Refs. 35 and 36). 

To study the various possibilities, a knowledge of the thermally- averaged 
dielectric properties of solid helium is essential. A calculation of the 
dielectric constant of solid helium, based on the self-consistent phonon 
formalism, is being undertaken for this reason. 

It is well known that electronic excitations, even the well localized 

ones, are able to migrate in crystals (see a discussion in Ref. 37). The 

mechanism for this migration is provided by resonant electronic interaction 

between neighboring excited and normal atoms. It is expected that in solid 

3 

helium the migration of 2 S excitations can well be represented by an effec- 
tive hopping Hamiltonian with parameters determined by evaluating 
thermally- averaged matrix elements of the He -He exchange interaction. 
Using this effective Hamiltonian, one should be able to compute the proba- 
bility of He -He collisions, for given He /He concentration, and thus to 
estimate the collisional loss rate of He atoms, on the assumption that a 

?! z 

near- neighbor collision of two He atoms is certain to lead to their deexcita- 
tion. In experiments with liquid helium, it has been found (Ref. 38) that 
activation of rotons and the presence of impurities independently tend to 
inhibit excited bubble collisions and so to lengthen the lifetimes of excited 
bubbles. It is interesting to speculate about possible ways of inhibiting 
"bubble 11 collisions in solid helium. 

V. METALLIC HYDROGEN 

In all cases where metallic hydrogen is believed to exist, i. e. , in the 
interior of Jupiter and in shock- wave experiments (Ref. 39), it is found at 
elevated temperatures, close to or above the melting temperature. Thus, 
from an experimental point of view, it is important to undertake studies to 
understand the physical properties of metallic hydrogen at finite, nonzero 
temperatures. These studies would also constitute a first step for future 
considerations of liquid metallic hydrogen. In view of the semiquantum 
nature of metallic hydrogen (Ref. 5), vacancies and interstitials are expected 
to have a dramatic effect on such properties of metallic hydrogen as lattice 
stability, melting, and superconductivity, to be discussed next. 
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A. LATTICE STABILITY AND MELTING 

The effect of static vacancies was considered by Aksenov (Ref. 9) who 
performed some simple model calculations for quantum solids. He obtained 
phonon softening due to vacancies. Also, in some of his calculations, lattice 
instability, characterized by at least one of the lattice frequencies becoming 
imaginary, occurred at vacancy concentrations as low as 5 to 10%. We 
intend to perform more realistic and complete calculations for metallic 
hydrogen, including also interstitials and their effect on lattice relaxation. 

It is known that in some cases (Ref. 40) interstitials cause a much larger 
lattice relaxation than vacancies, and therefore are of crucial importance 
in the study of defect effects. Not much is known about the relative impor- 
tance of vacancies and interstitials in metallic hydrogen, and this will be 
one of the first questions to be investigated in our studies. 

In many solids the lattice instability due to defects occurs at tempera- 
tures well above the melting temperature, and so is of no particular conse- 
quence. However, in metallic hydrogen there are indications that the 
situation is quite different. In a recent study of metallic hydrogen, Caron 
(Ref. 5) estimated a melting defect temperature, defined as the temperature 
at which the Boltzmann factor of the defect activation energy is of the order 
of 0 . 1 , corresponding to the maximum volume fraction of defects that the 
crystal can tolerate before vacancies and interstitials spontaneously 
recombine. For pressures of the order of the pressure at which a phase 
transition from molecular to metallic hydrogen is estimated to occur, the 
melting defect temperatures were found to be comparable to the ordinary 
melting temperatures of metallic hydrogen derived from other simpler 
models (Refs. 41 and 42). This implies that in metallic hydrogen the defects 
may well constitute the primary melting mechanism. 

Another aspect of finite-temper ature defects is worth noticing. While 
solid helium is a true quantum solid, metallic hydrogen is probably halfway 
between a quantum and a classical solid, as noted earlier. As a result, 
zero-temper ature defects are expected to be static, and phonons are prob- 
ably the only low-temper ature excitations. With increasing temperature 
both the number and mobility of defects increase, and it may be necessary 
to treat them as another type of excitation, similar to the vacancy waves 
(or vacancions) in solid helium (Refs. 10 through 13). 


JPL Technical Memorandum 33-701 


33 



B. SUPERCONDUCTIVITY 


Metallic hydrogen is believed to be a high-temperature superconductor 
(Refs. 3 and 4). The available estimates place the superconducting transi- 
tion temperature T c in the range 100 to 300 K. However, it should be 
noted that these estimates have been made on the basis of the simple BCS- 
type formula and McMillan's subsequent generalization (Ref. 43) for T c 
and they do not take into account either finite-temperature effects or the 
expected modifications caused by defect excitation. Now it is quite well 
established that defects and impurities in ordinary superconductors often 
have a drastic effect on T , either lowering it (usually) or raising it, and 
sometimes even causing oscillations in T c with increasing defect concentra- 
tion. In any reliable estimate of T c for metallic hydrogen, it is therefore 
essential to include the defect effects. A first guess might be that defects 
cause phonon softening in metallic hydrogen, thereby reducing the Debye 
temperature and the electron-phonon coupling constant, and thus, by the 
BCS formula for T , lead to a depressed transition temperature. With the 
detailed phonon spectrum to be evaluated, it will be possible to test this 
guess and, with the help of a theory of electron-phonon interactions for 
disordered structures (Ref. 44), to explicitly compute the changes of the 
superconducting transition temperature caused by the defects. 
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APPENDIX 


We provide here an integral representation for the variational thermo- 
dynamic potential, involving the full defect Hamiltonian {Eq. 37) rather than 
its linear approximation. Mathematically this involves replacement of the 
geometric series used in Eq. (39) by a sum of the form 



or, in real form 


2 


/ 


dt e 




1 - e a cos 2 Vb t 


1 + e 


-2a 


2e a cos 2 v^b” t 


(A- 2) 


For b - 0 this reduces, of course, to the geometric series. 

The thermodynamic potential is slightly more general, as it involves 
summation over the two variables N y = N 1 and N i = N 2 . The preceding pro 
cedure can be generalized to 
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a.N. - Y'' b..N.N. 

L~i i i L-t ij 1 j 


N 1 N 2 


N 1 N 2 


'f 


EiK 


dt„ e 


• t 2 - t 2 ^ {2a.t. + c.) N. - 2vt.t, 

1 2 / j ii l i Y l 2 


- K + t 2- ^1*2 


- I dt. dt- 

TT I 1 £ 


(l . e 2 " 11 * + C ‘) (l - e 2 “ 2 ‘ 2 + ° 2 


The five constants o\ , c., and v can be determined in terms of the five a., 

i l 1 i 

b_ by performing the t-integr ation in Eq. (A- 3a) and comparing with the 
original equations. The calculation is straightforward, though somewhat 
tedious, and gives 


c i = - a i • 


Y = - 


2 J h ll b 22 


2 _ b 12 ' 4b ll b 22 2 _ b l2 " 4b H b 22 

^ = 4b 22 ’ " 2 " 4b ll 


JPL, Technical Memorandum 33-701 


39 



Hence , 



(A-4) 


or, in real form. 



U 


A^Xj) (1 - A^x^) cosh 2yt^t^ - A. ^ i y 2 sinl1 Zytjtg 


(' - 2 Vl + A l) I 1 - 


2A 2 x 2 + A 2 J 


where 


- a. 



Finally, what is the difference between the full expression Eq. (A- 2) 
and the approximate value used in the text for b = 0? To answer this, we 
consider the difference between the two corresponding integrands 
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1 - e a cos 2 y/h t 1 

1 + e" 2a - e“ a cos zVht 1 - e“ a 

1 ^ (e~ a + e Za )(cos 2 \/b~ t - 1) _ f 

i -a - 2a ~ -a 0 [T~ . _ - a. 

1 - e 1 + e -2e cos 2 vb t 1 - e 

The approximation is justified if f « 1. For low temperatures, a (which is 
proportional to 1/kT by Eq. 39) is large and, consequently, 

r - a . - 2 a , 

f - e + e « 1 

For high temperatures, b is small (b is obtained from the quadratic terms 
in the full defect Hamiltonian Eq. (37)), cos (2 Vb t) - 1 ^ -2bt , and the 
contribution to the integral is again small. This brings us back to the low- 
defect approximation in which the b_ are neglected from the beginning. If 
necessary, however, one can start from the more rigorous Eq. (A-4). 
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